Constraining non standard recombination: A worked example 
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We fit the BOOMERANG, MAXIMA and COBE/DMR measurements of the cosmic microwave 
f~>) , background anisotropy in spatially flat cosmological models where departures from standard recom- 

■ bination of the primeval plasma are parametrized through a change in the fine structure constant a 

£\J \ compared to its present value. In addition to a we vary the baryon and dark matter densities, the 

spectral index of scalar fluctuations, and the Hubble constant. Within the class of models consid- 
ered, the lack of a prominent second acoustic peak in the measured spectrum can be accomodated 
either by a relatively large baryon density, by a tilt towards the red in the spectrum of density 
fluctuations, or by a delay in the time at which neutral hydrogen formed. The ratio between the 
second and first peak decreases by around 25% either if the baryon density flf,h 2 is increased or the 
, spectral index n decreased by a comparable amount, or if neutral hydrogen formed at a redshift z* 

^ • about 15% smaller than its standard value. We find that the present data is best fitted by a delay 

' in recombination, with a lower baryon density than the best fit if recombination is standard. Our 

best fit model has z, = 900, Q b h 2 = 0.024, Q, m h 2 = 0.14, H = 49 and n = 1.02. Compatible with 
present data at 95% confidence level 780 < z, < 1150, 0.018 < Q. b h 2 < 0.036, 0.07 < Q m h 2 < 0.3 
O ' and 0.9 < n< 1.2. 

o 
o 

Oh' 

I I. INTRODUCTION 

o 

Measurements of the Cosmic Microwave Background (CMB) anisotropy provide information about the physical 
conditions in the universe right before decoupling of matter and radiation, and allow for a precise determination of 
many cosmological parameters. The recently published results of the BOOMERANG JlJ and MAXIMA M experiments 
provide strong confirmation of the presence of the first acoustic peak in the anisotropy angular power spectrum at the 
degree scale, which suggests that the Universe is very close to spatial flatness. BOOMERANG and MAXIMA also 
provided high quality data at smaller angular scales, which do not reveal the presence of a prominent second acoustic 
- - 1 peak. 

The location, height and width of the first peak is in excellent concordance with the simplest spatially flat cosmolog- 
ical models with a nearly scale invariant spectrum of adiabatic density fluctuations (such as those motivated by generic 
inflationary models) and values of all cosmological parameters in good agreement with independent observations. 

The lack of a prominent second acoustic peak requires instead some sort of departure from the simplest models 
or most likely values of some cosmological parameters. Three possibilities have been highlighted i) the power 

spectrum of primordial density perturbations is tilted in the direction that suppresses small scale fluctuations; ii) the 
baryonic matter density is slightly above the upper value expected from big-bang nucleosynthesis; iii) the formation 
of atomic hydrogen was delayed by some mechanism that perturbed the standard ionization history of the Universe. 

The first mechanism reduces the height of the second peak simply because a spectrum tilted towards the red has 
less power at smaller angular scales. The second and third mechanisms work because the ratio between the heights 
of the second and first peaks decreases when the ratio R — jPb/p-y between baryon and radiation energy-densities at 
decoupling increases 0. Either if pb is increased or if recombination is shifted to lower redshifts R(z*) increases and 
the relative height of the second peak decreases. The shift of z* to lower redshifts also shifts the location of the peaks 
in the angular power spectrum towards larger angular scales. 

Peebles, Seager and Hu ^ developed a picture where very early sources of Ly a resonance radiation provide a delay 
in recombination rapid enough to avoid excessive dissipation of the first acoustic peak, and cause a 10% reduction in 
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the relative height of the second peak along with a 5% shift in the location of the first peak. Their model is essentially 
parametrized by the rate of excess of radiation by the early sources from that in the primeval plasma. 

In this paper we further explore the compatibility of present CMB anisotropy measurements with a non standard 
recombination history. We parametrize departures from standard recombination through a change in the fine struc- 
ture constant around decoupling compared to its present value. We then compute constraints in four cosmological 
parameters plus the fine structure constant from the COBE/DMR, BOOMERANG and MAXIMA data on CMB 
anisotropy. 

In section [il] we discuss the effects upon the CMB anisotropy spectrum of a variation in the fine structure constant 
at decoupling compared to its present value, and argue that they mimic the effects of other physical mechanisms 
that may change the recombination history. In section |l| we fit the COBE/DMR, BOOMERANG and MAXIMA 
data within spatially flat cosmological models allowing for variations of the fine structure constant, the baryon and 
dark matter energy densities, the Hubble constant, and the scalar fluctuations spectral index. We place bounds 
on the allowed variation of the fine structure constant, which translate into bounds on departures from standard 
recombination. In section [V we discuss the results. 



II. VARIATION OF THE FINE STRUCTURE CONSTANT AS A MODEL OF NON STANDARD 

RECOMBINATION 

Time-dependence in the fine structure constant a modifies the pattern of observed cosmic microwave background 
fluctuations Qualitatively, the main effect of a variation in a is due to the change in the binding energy of 

hydrogen. A change in a also changes the Thomson scattering cross section and modifies the recombination rates. 

We used the CMBFAST code with appropriate modifications to determine the effects of a change in a, along 
similar lines and with comparable results as in [B,[lO| . 
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FIG. 1. The visibility function in a model with standard recombination, and with recombination delayed (advanced) through 
a reduction (increase) of the fine structure constant a around decoupling by 3% 



The observed pattern of CMB anisotropies is largely determined by the visibility function, g(z) — exp{—T{z))d,T / dz, 
which measures the differential probability that a photon last scattered at redshift z. r is the optical depth due to 
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Thomson scattering. Of course, g(z) is extremely sensitive to the recombination history, since it largely depends on 
the time evolution of the fraction of free electrons. 

If a was smaller (larger) at recombination than its present value aojthe peak in the visibility function shifts towards 
smaller (larger) redshifts, and its width slightly increases (reduces) H[h|. These effects can be appreciated in Fig. [l], 
which displays the results of raising and lowering a by 3% around decoupling compared to its present value. 

The hydrogen binding energy scales asBoca 2 . Since the ionization fraction is largely determined by the Boltzmann 
factor exp (-B/T), the location of the peak of the visibility function should roughly scale as AT*/T* cx 2Aa/a in the 
limit of small Aa. Although this scaling is not exact ||, Az,/z, w 2Aa/a provides a reasonably good fit to the shift 
in the peaks in Fig. [j]. The width of the visibility function, 5z*, scales approximately as A(5z»)/<5z* w —Aa/a. 

Non standard recombination histories should have an effect upon the visibility function similar to that described 
above, except that the changes in the peak position and width may follow independent rates of change. In this respect, 
a variation of the fine structure constant alone can not mimic an arbitrary mechanism that may delay or advance 
recombination. At any rate, a variation of the fine structure constant is a good one-parameter emulator of realistic 
mechanisms for delayed recombination, since it incorporates the basic feature of shifting the decoupling redshift while 
maintaining a sufficiently fast recombination rate. Different physical mechanisms may delay recombination increasing 
the width of the last scattering surface at a different rate, compared to the shift in the peak, than a change in a. 
The important common feature is that the effect can be differentiated from changes in other cosmological parameters, 
such as the baryon matter density or the spectral index of density fluctuations. 

Indeed, the principal effect of a decrease in the value of the fine structure constant around decoupling is due to 
the shift of z* to lower values. The decrease of boosts the location of the acoustic peaks towards larger angular 
scales. A reduction of z* also implies an increase in i?(z*), the ratio between baryon and photon energy densities at 
decoupling, which increases the height difference between even and odd peaks. The height of the first peak is also 
affected by the change in the contribution of modes that entered the horizon during the radiation dominated period 
(the "early integrated Sachs- Wolfe effect"). Finally, the increase in the width of the visibility function that follows 
a decrease in the fine structure constant makes diffusion damping more effective, pushing down the damping tail of 
CMB anisotropics, while slightly increasing the still undetected degree of linear polarization [Q . This extra diffusion 
damping may help differentiate, when higher angular resolution data becomes available, a decrease of a or a delay in 
recombination from an increase in the baryon density or a red tilt in the spectral index. 

Time dependence of the fine structure constant is actually a theoretical possibility that deserves consideration by 
its own, not just as a model for delayed recombination. Indeed, unification schemes such as superstrings and Kaluza- 
Klein theories predict time variation of fundamental constants over cosmological timescales |l~3|-[ll|. More recently, 
a number of authors considered cosmological theories where the fine structure constant time dependence is due to 
the variation of the speed of light p6|-^8| (constraints on these theories from CMB anisotropy measurements were 
analysed in [fl9|| ). Furthermore, different versions of the above mentioned theories predict different time-dependence 
of fundamental constants. Thus, experimental bounds on their allowed variation are an important tool to check the 
validity of these theories. 

Constraints on the time variation of the fine structure constant have been placed from geophysical and astronomical 
methods. The Oklo natural nuclear reactor that operated about 1 .8 x 10 9 yrs ago in Oklo, Gabon ]2(J yields a constraint 
of —0.9 x 10~ 7 < Aa/a < 1.2 x 10~ 7 . Laboratory measurements based on comparisons of rates between clocks with 
different atomic numbers give a limit of Aa/a < 1.4 x 1CP 14 during 140 days pit]. From the analysis of natural 
long-lived a and (3 decayers in geological minerals and meteorites Dyson has placed a bound of Aa/a < 2 x 10~ 5 . 
The wealth of local tests, including possible correlated synchronous changes of different physical constants, lead to 
the estimate Aa/a < 2 x 10~ 5 for variations during the last few billion years p3|| . 

The astronomical methods are based mainly in the analysis of spectra from high-rcdshift quasar absorption systems 
p4] p7| . Most of the previous mentioned experimental data gave only upper bounds (e.g. the most stringent Aa/a = 
(—4.6 ± 5.7) x 10~ 5 for a set of redshifts z ~ 2 — 4 ^4|). Webb et al. ^5|, reported a positive measurement of the fine 
structure constant variation: Aa/a — (—1.09 ± 0.36) x 10~ 5 . 

Primordial nucleosynthesis provides a constraint to variations in a at the earliest times, derived form the relative 
abundance of 4 ife p8| . However, to compute this constraint a model for the a dependence of the proton to neutron 
mass ratio must be assumed. L. Bergstrom et al. (29J have arrived a the bound: Aa/a < 2 x 10~ 2 including in their 
analysis not only the A He abundances but also the abundances of other lighter elements that are much less model 
dependent. 
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III. DATA ANALYSIS 



We have performed a maximum likelihood analysis of the COBE/DMR, BOOMERANG and MAXIMA data 
within spatially-flat cosmological models with adiabatic density fluctuations. We computed the anisotropy cor- 
relation multipoles C; for a grid of models allowing for variations in a/ao = 0.86 — 1.04, Aa/ao = 0.02; 
uj b = 0.01 - 0.04, Au b = 0.003; u m = 0.05 - 0.6, Aui m = 0.05; Ho = 45 - 95, AH = 5 ; n = 0.8 - 1.15, An = 0.05. 
Here lj& = Vl b h 2 , uj m = Q m h 2 and ao is the present value of the fine structure constant. The cosmological constant 
was fixed to keep each model spatially flat, so that Oa = 1 — (o->b + Lo m )/h 2 . The Hubble constant Hq is measured 
in km/s Mpc, and h = H Q /(100 km/s Mpc). We included calibration errors (20% for BOOMERANG and 8% for 
MAXIMA) through the covariance matrix of measurement errors as described in |3(| . The amplitude A s of the scalar 
fluctuation spectrum was fixed minimizing \ 2 for each model. 

Our best fit model to the data, with 30 degrees of freedom reduced by 5 parameters, is characterized by: 

Lu b = 0.024, w m = 0.14, — = 0.91, H = 49 , n = 1.02 , x 2 = 30.55 (1) 

ao 

If we fix as a prior that recombination be standard (a = ao) the best fit becomes (with now 30 degrees of freedom 
reduced by 4 parameters): 

cj b = 0.031, w w = 0.19, H = 81, n = 1.02 , x 2 = 33.90 (2) 
If the prior is that the baryonic matter density ujb — 0.019 (the nucleosynthesis favoured value pl|) the best fit is: 

u) m = 0.1 , — = 0.9 , Ho = 45 , n = 0.96 , x 2 = 34.96 (3) 
a 

Fig. H displays the anisotropy angular power spectrum for the best fit model along with the COBE/DMR, 
BOOMERANG and MAXIMA data points, as well as the best fit model in the case of a standard recombination 
history. 
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FIG. 2. Best fit model allowing a variation of the fine structure constant (solid line) and best fit model with standard 
recombination (dashed line). 
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FIG. 3. 1-a and 2-a likelihood contours for Q.t,h 2 , fi m /i 2 , Ho and n vs. the relative change in the fine structure constant in 
spatially flat cosmological models. The delay in recombination is Az,/z* « 2(a/a — 1). 
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Two dimensional marginalized likelihood confidence contours of the four different cosmological parameters consid- 
ered as a function of the variation of the fine structure constant are shown in figure |. As we discussed in section 
H the lack of a prominent second acoustic peak in the BOOMERANG and MAXIMA data can be accommodated 
within simple spatially flat cosmological models either with a relatively large baryon density, a tilt towards the red in 
the power spectrum, or a delay in recombination. Our contours show that current data favour delayed recombination 
over the other options, which at any rate are all consistent at the 95% confidence level. 

The likelihood contour in the plane (cob, a) reveals a degeneracy between reducing the baryonic matter density 
and delaying recombination. The best fit model, in which a is smaller by more than 8% than its present value, and 
thus recombination is delayed by around 16% in redshift, has indeed a baryon density more than 25% smaller than 
the best fit model with standard recombination. Still, the best fit model has a baryon density slightly above the 
nucleosynthesis upper bound. 

The best fit when the baryon density is fixed to the nucleosynthesis preferred value (wf, = 0.019) requires a large 
(20%) delay in the recombination redshift, in addition to a small (about 4%) tilt of the spectral index towards the 
red. 

The two-dimensional contours also reveal that fit of the data in delayed recombination scenarios is best with 
relatively small values of the Hubble constant. The degeneracy with the value of the spectral index is instead rather 
small. All the range allowed for departures from standard recombination at 95% confidence level is compatible with 
a scale invariant (n — 1) spectrum. 
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FIG. 4. One dimensional marginalized likelihoods. Left panel: marginalized over u>b, u) m , Ho and n, as a function of a/ao 
(without priors) and with uii, fixed to its nucleosynthesis prediction; right panel: marginalized over Lu m ,Ho, n and a as a 
function of uib (without priors) and with a fixed to its present value ao (standard recombination) . 

One dimensional marginalized likelihoods are shown in figure |]. These make clearer that, within the family of 
cosmological model considered, the present data are best fit in a delayed recombination scenario. If the baryon 
density is within the nucleosynthesis bound then delayed recombination is imperative (a tilt in the spectral index 
alone is not enough). 

At 95 % confidence level, the bounds on the constrained parameters are 

0.018 < n b h 2 < 0.036 , 0.07 < n m h 2 < 0.3 , 0.86 < a/a Q < 1.03 , H < 95 , 0.9 < n < 1.2 (4) 
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We can use the CMB anisotropy data to place a bound, within the family of cosmological models considered, upon 
departures from standard recombination. We computed z* (determined as the location of the peak in the visibility 
function) for each model in the grid. Then we found the model with minimum \ 2 f° r z * m given fixed intervals. From 
this result we built a one-dimensional likelihood for z* , shown in figure ^[ 
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FIG. 5. One dimensional marginalized likelihood for z* 
From figure || we estimate the allowed range for z* to be, at the 95% confidence level, 

780 < z» < 1150 (5) 

with a best fit value around z* = 900. 

Notice that the result in eq. || is in very good agreement with the range in z» one would expect from the range 
0.86 < a/ao < 1.03 and the approximate scaling z* = 1080[1 + 2(a — ao)/«o]- 

IV. CONCLUSIONS 

We have modeled non standard recombination through a variation of the fine structure constant a around decoupling 
compared to its present value ao- We performed a maximum likelihood analysis of the BOOMERANG, MAXIMA and 
COBE/DMR measurements of the cosmic microwave background anisotropy in a grid of spatially flat cosmological 
models with five parameters varied independently: a, Sl;,/i 2 Q m h 2 , h and n. The recombination redshift scales 
approximately as z* = 1080[1 + 2(a — ao)/ao\. 

The lack of a prominent second acoustic peak in the data suggests that the best fit should have, within the family 
of models considered, either a relatively large baryon density compared to the nucleosynthesis value cut, = 0.019, a 
tilt away from scale invariance towards the red in the spectrum of density fluctuations, or a delay in the time at 
which the CMB decoupled from matter. The best fit is that recombination was largely delayed to a redshift around 
z* = 900, with a nearly scale invariant spectrum of fluctuations, a low value of the Hubble constant, H = 49, and 
a baryon density f^/i 2 = 0.024. A delay in recombination allows a good fit to the data with a lower baryon density 
as compared to a standard ionization history. Acceptable fits with a baryon density as low as il^h 2 = 0.018 can be 
achieved with longer recombination delays. Good fits with such low baryon densities can not be achieved with just a 
tilt in the spectral index and a standard recombination history. Recombination can not be delayed beyond z* w 780 
at 95% confidence level and within the cosmological models considered. Notice that in our grid of models Ho > 45. 
Recombination may eventually be delayed to lower values of z* at the same confidence level if this (realistic) prior 
assumption upon Hq were lifted. 

It is instructive to discuss these results from a semianalytic approach. It has been shown f32|| that the power 
spectrum of CMB anisotropies up to the third acoustic peak can be conveniently characterized by four observables, 
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namely: the position of the first acoustic peak li, the height of the first peak relative to COBE normalization Hi, the 
height of the second peak relative to the first H2 and the height of the third peak relative to the first H3. We have 
performed a semianalytic fit around our best fit model of eq. ([!]) assuming a linear dependence on the parameters, 
with the approximate result: 

^1 = 1.4^ + 0.05^ - 0.14^ - 0.17^ + 0.19An (6) 

h ol uj b uj m Hq 

= 0.9^ + 0.44^ - 0.47^ - 0.13^° + 2.6An (7) 

Hi a LU b 0J m H Q 

^1 = 3.2^ _ 0.85^ - 0.04^ + 0.96An (8) 

H 2 a u b uj m 

^ll = 2 .3^ - 0.35^ + 0.38^ + 1.2An (9) 

H 3 a u b uj m 

The peaks location and heights of the best fit model are given by l\ = 206.5, H% — 7.8, H% — 0.37, H3 = 0.43. 

Equation (|J) condenses most of the discussion we already made in sections | and || about the effects that can 
accommodate the absence of a prominent second peak in the present CMB anisotropy data. Indeed, the relative height 
of the second acoustic peak compared to the first decreases by about 10% if a decreases around 3%. Alternatively, 
the same reduction in H2 can be obtained by an increase of around 10% in the baryon density, or a 10% tilt towards 
the red in the spectral index of density fluctuations. H2 is instead quite insensitive to u> m and Ho- 

Notice from eq. (jfy that a decrease in a by 3% not only reduces H2 by around 10% but also shifts the location of 
the first acoustic peak by almost 5% towards larger angular scales. This coincides with the results of Ref. ||, for a 
specific model of delayed recombination based in early Ly a sources of ionizing radiation. 

The data from BOOMERANG and MAXIMA does not have the accuracy and angular resolution to test for the 
presence of a third acoustic peak. Consequently, we should not use H% as a relevant parameter. Eqs. (^-|^) alone 
are insufficient to analyse degeneracies as a function of just one parameter. It is however quite apparent from the 
two-dimensional contours of the previous section that the weakest degeneracy is that on the spectral index n, which 
is quite close to scale invariance. If we assume An = we derive, imposing AZi = AHi = AH2 = 0, 

A^ = 36 Aa ; A^ =35 Aa j ^ = M Aa 
Ub a uj m a Ho a 

which is in very good agreement with the shape of the two-dimensional likelihood contours of the previous section. 

Variations in u> b , a and n modify i?2- Neither the variation of u> b nor of n have a strong impact upon the location 
of the first acoustic peak. Variation of a does, as we already discussed. Compensation of the peak shift as a decreases 
requires a large decrease either of H or of oj m . Something similar happens with Hi , which is also very sensitive to 
changes in the spectral index n. Our results reveal a strong degeneracy in the plane (ck,Hq), which appears to be 
quite robust, in the sense that it persists with similar strength even under variations of n. The degeneracy in the 
plane (a, ui m ) is instead quite sensitive to the precise degeneracy in the (a, n) plane. Information on the structure of 
the angular power spectrum of CMB anisotropies around the location of the third acoustic peak is necessary to more 
accurately pinpoint this issue and eventually break some of the present degeneracies. 

A decrease in a, as well as any delayed recombination mechanism in which there is an increase in the width of 
the visibility function, makes diffusion damping more effective, which increasingly reduces the heigths of the peaks at 
smaller angular scales. This is already noticeable around the third acoustic peak. A decrease in a that reduces H2 as 
efficiently as an increase in tu b has a stronger effect in reducing H3, as can be seen from eqs. (^,^|). Notice, for instance, 
that the best fit model with delayed recombination of Fig. g has a third accoustic peak suppressed compared to the 
best fit with standard recombination, which has a larger baryon density. Measurements of the spectrum around and 
beyond the third accoustic peak can potentially break the degeneracy between the effects of an increase in uj b and 
delayed recombination. 

Our results can also be seen as a bound on departures from standard recombination, which can happen at redshifts 
as low as z* = 780 at 95% confidence level, or as a bound on the time variation of the fine structure constant between 
decoupling and the present time. The bound we obtain, 0.02 > Aa/a > —0.14, is an order of magnitude weaker 
than expected in j| from future satellite measurements. The weakness of this bound is a consequence of the large 
present degeneracies. More accurate measurements and higher angular resolutions will certainly improve this bound. 
Whether a full order of magnitude improvement is or is not achieved will depend on how much this degeneracy is 
actually broken. 

Our result that there is a preference for a smaller value of a at decoupling compared to its present value is in 
agreement with similar conclusions of recent work p^,p3] that also analysed BOOMERANG and MAXIMA data 
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allowing for a variation in the fine structure constant. In the present work we used the time varation of the fine 
structure constant as an example of how to model non-standard recombination histories. A preprint that develops 
a two-parameter model for non-standard recombination histories |55| appeared after submission of our paper, with 
compatible results. 
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